Creation and detection of skyrmions in a Bose-Einstein condensate 
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We present the first experimental realization and characterization of two-dimensional skyrmions 
and half- skyrmions in a spin-2 Bose-Einstein condensate. The continuous rotation of the local spin 
of the skyrmion through an angle of 7r (and half-skyrmion through an angle of 7r/2) across the cloud 
is confirmed by the spatial distribution of the three spin states as parameterized by the bending 
angle of the ^-vector. The winding number, w = (0,1,2), of the skyrmions is confirmed through 
matter- wave interference. 
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Skyrmions are topological solitons that were first envi- 
sioned in the 1960's as part of a nonlinear field theory to 
model mesons and baryons in nuclear physics [1]. They 
are one of the many topological defects that can be de- 
scribed and classified by homotopy groups of their order- 
parameter space 0, El- A 3D skyrmion is an S 3 — » S 3 
map that can be labeled by an integer representing the 
topological degree, or winding number, of the map [![• 
The 2D or "baby skyrmions" are characterized by a local 
spin that continuously rotates through an angle of tt from 
the center to the boundary of the system [5[ , and are the 
focus of this Letter. The search for skyrmions spans sev- 
eral fields including superfluids, solid state physics, liq- 
uid crystals, and superconductors 

The experimental results presented in this Letter mark 
the first deterministic creation of skyrmions and half- 
skyrmions in a spin-2 Bose-Einstein condensate (BEC). 

As in many condensed matter systems, a BEC of al- 
kali metal atoms can be described by an order parameter 
which becomes vectorial with 2F + 1 components when 
multiple states in the same hyperfine spin-F manifold are 
simultaneously populated [13[ . These spinor condensates 
have made it possible to study experimentally a multi- 
tude of spin excitations such as coreless vortices, solitons, 
and spin domains in a dilute system instead of one dom- 
inated by the strong interactions and high densities of 
condensed matter physics 
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For a spin-2 condensate, two five- component order pa- 
rameters are required to sufficiently define the system: 
(/), which describes the ferromagnetic character, and 6, 
which is related to the dynamic creation of spin singlets 
through collisions [HI, [ljj]. In our experiments, prior to 
the creation of the skyrmion, the BEC is spin-polarized in 
the \F = 2, m F = 2) = |2) state with a density of « 10 11 
cm -3 . This means that initially (/) = 2, while 6 is neg- 
ligible for the timescales involved. 

A detailed description of our coherent Raman interac- 
tion can be found in [2(j. At the start of the experiment, 
the BEC is untrapped and cylindrically symmetric with 
a magnetic field of B = 1.33 G oriented along the z-axis. 



Two beams (<j~ , <r + polarized) propagate collinearly and 
parallel to the quantization axis, creating a double A (or 
"M") system by simultaneously coupling the |2) state 
with the |2,0) = |0) and |2, -2) = |-2) ground states. 
The Raman beams are applied diabatically to the BEC 
in 5 /is pulses. The short timescale enables the creation 
of non-equilibrium spin textures in the condensate, while 
the co-propagating beam geometry minimizes the trans- 
fer of linear momentum. 

The cr~, <j + Raman beams have first-order Laguerre- 
Gaussian (LG -1 ) and Gaussian intensity profiles, respec- 
tively, so the population transferred to the |0) (|— 2)) 
state acquires a w = 1 (w = 2) azimuthal phase winding. 
A core of w = atoms is left in the initial |2) state cor- 
responding to the intensity minimum of the LG beam. 
The relative population transfer to the |0), |— 2) states 
is a function of the two-photon detuning of the Raman 
beams [20|, while the spatial dependence of the three 
components of the coreless vortex additionally depends 
on the beamwaist of the LG beam. It is in this way that 
the Raman beam coupling, combined with the relative 
energetic degeneracy of the |2) <-> |0) and |0) <-> |— 2) 
transitions, creates a coreless vortex with winding num- 
ber: w = (0, x, 1, x, 2). The x indicates that there is no 
population in a spin state, and so there can be no phase 
winding there (the x in the |2,1), |2, — 1) states will be 
dropped) . The Raman interaction effectively evolves the 
order parameter of the initially spin-polarized BEC to 
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Here n(r) is the density of the cloud, (r, ip) are polar 
coordinates, and /3(r) is the bending angle which charac- 
terizes the rotation or "bending" of the local spin across 
the cloud. The local spin is represented by the spin tex- 
ture's ^-vector: 

£(r, (p) = z cos (3(r) + sm(3(r)(x cos cp + y sirup). (2) 



This ^-vector description of the local spin is a general 
definition of a 2D skyrmion, and it requires that the 
bending angle be monotonic and satisfy /3(0) = and 
/3(R) — 7r, where R is the boundary of the cloud [5[. 
For the half-skyrmion, the ^-vector sweeps through half 
the angle of the skyrmion, and so 3 must monotonically 
satisfy /?(0) = 0, (3(R) = tt/2 @,[2lj]. The continuous t- 
vector and non-singular order parameter are fundamen- 
tal properties of skyrmions and will be discussed in more 
depth later. 

Previous experiments conducted in spin-2 87 Rb have 
focused on equilibrium studies and spin-mixing dynam- 
ics without vortices [22[ . Initial coreless vortex work has 
involved coupling the F = 1, 2 ground state manifolds of 
87 Rb by RF transitions, and the adiabatic manipulation 
of spin-1 23 Na by magnetic fields 3, 15]. While novel in 
their approach, these coreless- vortex techniques are lim- 
ited in their applicability: the former by its inability to 
create mult i- component coreless vortices within a single 
ground state manifold, and the latter due to the inherent 
difficulty of creating and controlling arbitrary magnetic 
fields. In contrast, we can tailor the Raman interaction 
used here to create w = (0,1,2) coreless vortices while 
simultaneously controlling the spatial distributions of the 
spin states and the normalized magnetization of the cloud 
in order to precisely engineer complicated spin textures 
such as skyrmions (Fig. [T]) and half-skyrmions (Fig. [2]). 

The (0, 1, 2) coreless vortex presented in Fig. [T] is a 2D 
skyrmion created in a spin-2 BEC. The absorption im- 
age, DJa), shows the three components of the cloud after 
they have been spatially separated by an inhomogeneous 
magnetic field for imaging. Prior to this "Stern- Gerlach" 
pulse however, the coreless vortex is axi-symmetric and 
the cloud maintains its cylindrical symmetry (Fig. 0Jc)). 
The spatial dependence of the spin states comprising the 
spin texture, Fig. QJb), can be used to find the bending 
angle of the cloud by fitting (pQ) to lineouts taken from 
the absorption image. The solid lines in DJb) are the re- 
sult of simultaneously fitting all three spin states to $1]) 
for a single bending angle, f3\. The agreement confirms 
that this spin texture is indeed a skyrmion. 

The polarization of the cloud, Fig. QJd), is one of 
the most recognizable differences between skyrmions and 
half-skyrmions. If the density of the cloud were constant, 
the polarization, or local magnetization, of the skyrmion 
would be equivalent to the z-component of the i?- vector: 
i z = cos/3 [2l(. Since this is not the case for a BEC, 
the plot of t z has been multiplied by the density profile 
of the cloud. This necessarily supresses the amplitude 
of the polarization at the boundary, but in the outer re- 
gion of the cloud the dominance of the |— 2) state still 
causes the polarization to become negative. It should 
be noted that the solid line is not a fit to the polariza- 
tion data, instead it is a curve created using the bending 
angle generated by fitting ([T|) to the absorption image. 
The close agreement confirms that f3\ parameterizes the 
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FIG. 1: The absorption image (a) of a 2D skyrmion created 
in spin-2 87 Rb. The winding number, w, for each spin state 
is indicated. The contrast of the image has been increased to 
make the |— 2) state more visible. The relative populations 
of the |2), |0), and |-2) states are: 51%, 30%, 19%. (b) Az- 
imuthally averaged lineouts (points) for each spin state agree 
well with the plots of (pQ) (solid lines) for the f3± (f) determined 
by a nonlinear least-squares regression, (c) 3D plot of the solid 
lines in (b), where the colors red, green, and blue correspond 
to the 1 2), |0), and | — 2) states, respectively. The winding 
number of the spin state is indicated by the number of ar- 
rowheads, (d) The polarization of the skyrmion. The points 
represent \4>2(r)\ 2 — |0_2(f)| 2 from (b), and closely match the 
solid line: i z (r) = n(r) cos /3i (f). (e) The vector plot of the 
skyrmion ^-vector clearly shows the rotation of the local spin 
through an angle of tt across the cloud. The same (3\ was 
used for all theory plots, and the goodness of fit for (b), (d) 
is Rj = 0.977, Rj = 0.982, respectively. 



cloud well, and can therefore be used to reveal the spatial 
dependence of the local spin of the cloud. 

Fig. QJe) presents the ^-vector of the skyrmion: ini- 
tially parallel to the z-axis at the center of the cloud, it 
continuously rotates through an angle of tt to lie anti- 
parallel to the z-axis a,t f = r/R = 1. This is a cylin- 
drically symmetric 2D spin texture, analogous to the 
Anderson- Toulouse spin texture predicted to exist in su- 
perfluid helium In contrast to the skyrmion ground 
states predicted to arise spontaneously in magnetic met- 
als [sl , the rotation of the i?- vector is in the r — z plane. 
This spin texture is reminiscent of the spin helices re- 
cently observed in semiconductor quantum wells [8[ , and 
those predicted to occur in nematic liquid crystals with 
phase angle a = [10]. Creating a skyrmion in a BEC re- 
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quires a coreless vortex where the successive peaks in the 
spatial distributions of the components contribute to the 
continuous bending of the ^-vector, while simultaneously 
ensuring that the order parameter is non-singular. 
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FIG. 2: Absorption image (a) of a half-skyrmion created in 
spin-2 87 Rb. The lineouts (b), reconstruction (c), polarization 
(d) and ^-vector (e) have been fit, processed and presented as 
in Fig. [1] The relative populations of the |2), |0), and |— 2) 
states are: 54%, 31%, 15%, and the goodness of fit for (b), (d) 
is R\ = 0.996, Rl = 0.999, respectively. 



The half-skyrmion or meron spin texture was first en- 
visioned in the 1970's, and was predicted to form sponta- 
neously in supernuid helium [7U2lLl23||. Fig. [2] shows the 
first confirmed creation of a half-skyrmion. This 2D spin 
texture can also be described by (pQ), but the bending 
angle must monotonically approach tt/2 at the boundary 
of the cloud instead of tt. This has a profound effect on 
the spatial distribution of the three components of the 
coreless vortex, as can be seen in[2fb,c) and the resulting 
spin texture, (e). 

The polarization, [2jd), of the half-skyrmion is quali- 
tatively distinct from that of the skyrmion. Maximum 
at the origin, it monotonically tends to zero across the 
cloud. This requires \4>2(r)\ 2 and |</>_2(r)| 2 to decay to 
zero together as r — > R, and places an upper bound on 
the relative population of the |— 2) state. These features 
are clearly present in [2][b) and particularly (c), where 
the presence of all three spin states (represented by red, 
green, and blue) at the boundary of the cloud results in 
a turquoise color instead of the solid blue (indicative of 
the |— 2) state) of the skyrmion (Fig. Hfc)). 



By fitting (pp) to the absorption image lineout data in 
Fig. [2](b), we again find a bending angle that parameter- 
izes the spatial dependence of the spin states. The result- 
ing ^-vector of the half-skyrmion is shown in Fig. [2je). 
The local spin of the cloud still points along z at f = 0, 
but in contrast to the skyrmion, rotates only through an 
angle of tt/2 to lie in the x — y plane at the boundary. 
The effect of the monotonicity constraint on the bending 
angle parameter is more obvious for the half-skyrmion. 
At any point where (3 > tt/2 the local spin would be neg- 
ative, violating the definition of a half-skyrmion. Recent 
theoretical work in high temperature superconductivity 
has featured the half-skyrmion [l2j], but it has yet to be 
observed experimentally in those systems. 
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FIG. 3: Matter-wave interference of (0,1,2) and (0,-1,-2) 
coreless vortices alongside numerical predictions. The con- 
trast of the data and corresponding theory has been increased 
for the |0), | — 2) states to make the structure more visible. The 
relative populations of the |2), |0), |— 2) states for (a) and (b) 
are: 49%, 30%, 21% and 50%, 30%, 20%. 



The azimuthal phase winding of the spin states com- 
prising the skyrmions is critical to their topological sta- 
bility. We use matter-wave interference to directly con- 
firm that the Raman interaction is creating the (0, 1,2) 
coreless vortex structure. The diabatic pulses and non- 
equilibrium nature of the interaction make it possible to 
generate coherent superpositions of coreless vortex states 
through repeated application of the Raman beams [241 ]. 
By changing the order of the a~ beam from LG _1 to 
LG +1 between the pulses, we create a superposition of 
a (0,1,2) and (0,-1,-2) coreless vortex. The interfer- 
ence patterns revealed upon absorption imaging confirm 
the vorticity of the spin states [25[. Fig. [3] (a) and (b) 
show absorption images of such superpositions, created 



4 



under experimental conditions simliar to the skyrmion 
and half-skyrmion presented in Figs. [T] and O respec- 
tively. The two azimuthal nodes in the |0) state confirm 
that the interfering clouds have orthogonal azimuthal 
phase windings: 1 — (—1) = 2, while the clover leaf in- 
terference pattern in the |— 2) state is clear evidence of 
the superposition of a w = +2, w = — 2 vortex state. 
This simultaneously demonstrates the stability of these 
spin textures over the 20 ms of time-of-flight expansion 
between their creation and imaging. 

Skyrmions can be created across a range of magneti- 
zations. The normalized magnetization is critical in de- 
termining the stability and subsequent evolution of spin 
textures [JJJ |2l|, [26[ , as it is conserved for many interac- 
tions f22(| . For a spin-2 system, M/N ranges from [—2, 2], 
and the skyrmions presented in Figs. HJ El have normal- 
ized magnetizations of M/N — 0.64 and 0.78, respec- 
tively. It is unknown at present if the ground state phase 
in zero magnetic field for spin-2 87 Rb is polar or cyclic 
[Hj]. However, the range of normalized magnetizations 
for a (0,1, 2) coreless vortex to be stable in polar spin-2 
87 Rb is estimated to be [0.25, 1.3] [19], and so from this 
perspective we would expect both of these spin textures 
to be stable. The implementation of an optical dipole 
trap is underway to enable studies of the stability and 
evolution of not only the skyrmions presented here, but 
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also other non-trivial coreless vortex states 

In this Letter we have presented, in the form of (0, 1, 2) 
coreless vortices, the first experimental realization of 
skyrmions and half-skyrmions in a spin-2 condensate. 
The coreless vortices were shown to be well-characterized 
by bending angles, ^1,2 (r), that enabled the reconstruc- 
tion of the ^-vectors of the respective 2D spin textures. 
This revealed not only the continuous evolution of the 
local spin of the clouds, but also their non-singular pro- 
files. Further investigation into the stability of skyrmions 
in spinor condensates has the potential to strengthen the 
connections between BEC physics and other systems with 
spin degrees of freedom [28[ . 
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